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The possible effect of environment on the efficiency of a 
quantum algorithm is considered explicitely. It is illustrated 
through the example of Shor's prime factorization algorithm 
that this effect may be disastrous. The influence of environ- 
ment on quantum computation is probed on the basis of its 
analogy to the problem of wave function collapse in quantum 
measurement. Techniques from the Hepp-Colemen approach 
and its generalization are used to deal with decoherence prob- 
lems in quantum computation including dynamic mechanism 
of decoherence , quantum error avoiding tricks and calculation 
of decoherence time. 

PACS number(s): 03.65-BZ, 71.10.Ad,89.70.+c 

1. Introduction 

Quantum computations (QC) can be understood as 
a quantum- mechanical time evolution of certain quan- 
tum systems (so-called qubits) [1-9], in which the non- 
classical dynamic feature, such as the quantum coherence 
of states, plays a dominant role. Indeed, it is the purely 
quantum characters that makes it possible for a theo- 
retical quantum computer to solve certain hard math- 
ematical problems efficiently. In this respect, perhaps 
the most important example is Shor's prime factoriza- 
tion algorithm [6]. As the quantum computation is a 
quantum process, preserving coherence, at least to some 
extent, throughout the whole process is thus an essen- 
tial requirement. In fact, the decoherence resulting from 
the coupling with environment may make an quantum 
algorithm invalid and may cause unwelcome exponential 
increase of errors in output results [10-12]. Actually, the 
decoherence process was even regarded as a mechanism 
for enforcing classical behavior in the macroscopic realm 
[13]. In this view a decohered quantum computer be- 
comes a classical one. To overcome this difficulty caused 
by decoherence, some schemes have been proposed in the 
last several years [14-20]. Among them are the quantum 
error correcting technique inspired by the classical error 
correction theory and the schemes avoiding decoherence 
presented in connection with the strategy preventing de- 
coherence from certain subsets of quantum states [21-23]. 
The latter was systematically described in the framework 
of error avoiding quantum coding (EAQC) . For the math- 
ematical details we refer the readers to ref.[23]. 

It is recongnized that there exists a substantially close 



relation between the problem of decoherence in quan- 
tum computation and the problem of wave function col- 
lapse(WFC, also called von Neumann's reduction) in 
quantum measurement [24]. In the view point of quan- 
tum dynamics [25-35] for quantum measurement, both 
the measured system and the measuring instrument (de- 
tector) obey the Schrodinger equation and the dynamic 
evolution governed by their interaction can result in the 
WFC under certain conditions. For example, the collapse 
happens if the detector contains a great number of par- 
ticles or if the detector is in a state with a very large 
quantum number. These two cases are usually referred 
to as macroscopic limit and classical limit respectively 
[30] . We recall that in the traditional theory of quantum 
measurement [24], the WFC postulate is only an extra 
assumption added to the ordinary quantum mechanics. 
Under this postulate, once we measure an observable and 
obtain a definite value Ofc the state of the system must 
collapse into the corresponding eigenstate |fc) from a co- 
herent superposition \<t') =^f;Ck\k){k\. In the terminol- 
ogy of density matrix this process can be described by a 
projection p = \4>){4>\ p = from a pure 

state to a mixed state. This projection process means the 
loss of quantum coherence. There is a strong resemblance 
between this process and the quantum decoherence of a 
quantum computer resulting from the coupling with the 
surrounding environment. Thus, with the correspondence 
of the environment surrounding the quantum computer 
to the measurement instrument monitoring the measured 
system, some obtained results in the quantum dynamic 
models [26-36,] based on Hepp-Colemen (HC) approach 
[25] for quantum measurement can be applied to discuss 
decoherence problems in quantum computation such as 
the dynamic mechanism of decoherence in quantum com- 
puter, the quantum error avoiding techniques and the 
calculation of the decoherence time of a quantum com- 
putation process. We can also re-consider the strategy 
grouping the quantum states of qubits to form decoher- 
ence free subsets [21-23] and analyse the dynamic process 
decohering the qubits beyond the decoherence free sub- 
sets . 

In comparison with the dynamic theory of quantum 
measurement, the system of qubits in quantum compu- 
tation is an open system S (similar to the measured 
system) surrounded by an environment ( similar to the 
measuring instrument, detector). The environment may 
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cause two unwelcome effects on computation process, 
namely, causing the states of the qubit system to con- 
tinuously decohere to approach classical states [13, 25- 
36] and dissipating the energy of the qubit system into 
the environment. Mathematically, they are respectively 
described by vanishing of the off-diagonal and diagonal 
elements of the reduced density matrix of 5* . The dissi- 
pation effect of imperfect isolation happens at the relax- 
ation time-scale Trei ■ It is relatively easy to make systems 
having a very large Trei and thus allowing a reasonable 
number of operations to complete [10]. In contrast, the 
effect of decoherence, is much more insidious [10, 37-42] 
because the coherence information leaks out into the en- 
vironment in a time scale much shorter than 
quantum system evolves [10, 37-42]. In fact, it is rather 
difficult to realize a quantum system with the time scale 
Td for decoherence smaller than that for the dissipation. 
For example, let us consider an oscillator in a superposi- 
tion of two coherent states separated by a distance I from 
each other. Suppose they interact linearly with a bath of 
oscillators at temperature T. Then the decoherence time 
for this system is linearly proportional to the relaxation 
time with a rather small ratio ( )^ at high tempera- 
tures where Xthem is the thermal de Broglie wavelength 
[10, 37-42]. This implies that the qubits decohere much 
faster than they dissipate. Thus, the sensibility of quan- 
tum computation mainly depends on r^rather than Trei- 
For this reason , the present discussions in this paper 
only focus on the decoherence problem. 

The arrangement of this paper is as follows. In sec. 2, 
we consider the influence of environment on Shor's prime 
factoring algorithm; From sec. 3 to sec. 7 we consider de- 
coherence problem caused by environment with concrete 
models; In sec. 9, we discuss the universality of environ- 
ment in the weakly coupling limit; Finally in sec. 8 we 
draw some conclusions. 

2. Decoherence in Shor Factoring Algorithm 

In this section we illustrate the possible influence of 
environment on the validity of a quantum algorithm 
through an example, Shor's prime factorization algorithm 
. In this case, we recall, the so called quantum computer 
has two registers in use. According to Shor's method [6] , 
to factor a number n one should flrst of all choose a num- 
ber X . Then the first step is to put the flrst register in 
the uniform superposition of q states \a){a = 0, 1, ....g— 1) 
and the second one in a single state |0). This leaves the 
machine in the state 

\m) = ^J2\a)^\0) (2.1) 

Next, one computes x'^mod.{n) in the second register, 
leaving the machine in the state 
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= — ^ |a, x''mod.{n)) (2.2) 
v9o=o 

Then one performs a Fourier transform Aq on the first 
register. This leaves the machine in the state 

= i|]|]exp[^]|c,a;«mod.(n)) (2.3) 

^ 0=0 c=0 ^ 

Finally one observes the machine. One easily finds that 
the probability that the machine ends in a particular 
state \c,x^mod.{n)) = \c,k{x;n)) is 

P(c,fc) = ^| Yl exp[^]|2 (2.4) 

a:x'^=x^Tnod.{n) 

Shor shows that if c lies in a particular region one can 
determine a nontrival factor of n from the value of c. 
Denote the one-try-success probability of this method by 
Ps- Then one has the following result: 

Ps > r(t>{r)p{c, fc) > 1/ In n > 1/3 In n (2.5) 

where r is the least integer such that x^ = l{mod.n) and 
(p is Euler's totient function. 

Now let us take the influence of environment into ac- 
count to some extent. Asume that the environment 
is comprised by N particles. In this case we denote 
by (t)'{0),(j)'p{t),p'{c,k) and p'^ the correspondences of 
<f){Q),(f)F{t),p{c,k) and ps respectively. Then we have 

V« a=0 

where |e) = jei) (g) [62) ® ... ® jeAr) is the initial state of 
the environment without correlation with the state of the 
machine. Here, \ek){k = 1,2,.. .,N) denotes the initial 
states of individual particle comprising the environment 
. Accordingly, 

1 

\(l)'{t)) = — ^ |a) O Ix^mod.in)) \e[a]) 
v9 0=0 

1 

= — V |a, x°-mod.{n)) \e[a]) (2.6) 

V^a=0 

where |e[a]) = Ua{t)\e) and Ua{t) is the effective evo- 
lution operator of the environment correlated with the 
state I a). For simplicity we do not consider the influence 
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of environment in the process of the Fourier transform 
^g.Thus we have 

\<kW)) = - I] I]exp[?^] |c, x''mod.{n)) ® \e[a]) 

1 a=0 c=0 ^ 

As the only difference between the present model and the 
original one is the involvment of the environment vari- 
ables \e[a]) in the entanglement, to proceed along with 
the discussion we should consider the reduced density 
matrix 



p{t)=Tr,{\Mt)){<t>F{t)\) = 



quantum measurement as a consequence of a certain fac- 
torizable structure of the effective evolution operator of 
environment [30-36]. In fact, if Ua can be factorized as 



(2.12) 



where Ul{t) only concerns the j'th particle in the envi- 
ronment, the decohering factor can be expressed as A''- 
multiple product 



N 



F{a,a')=J{{e^\Ui){t)Ui{t)\e^) 



^ 9-1 9-1 (J-l 9-1 
^ 0=0 c=0 a'=0 c'=0 



2'jTi{ac — a'c'). 



N 



^l[\F^{a,a') 



(2.13) 



{e[a']\e[a])\c,x°'mod.{n)){c',x'''mod.{n)\ (2.8) 

Here wc have traced over the environment variables. We 
notice that the contribution of environment is given by 
the transition matrix element 



F{a,a') ^ {e[a']\e[a]) = (e|C/t,(t)f/„(i)|e) (2.9) 

F{a,a') is usually called decohering factor. Now it di- 
rectly follows that 

p'{c,k)=Tr{p{t)\c,k){c,k\} 



9-1 

E 



g2 



2m{a- a')c , 
exp[ —]F{a, a ) 



(2.10) 



We are now in a position to consider two extreme cases. 
For the first case, suppose that the qubit system is com- 
pletely isolated. In this case we have Ua' = Ua for a' ^ a, 
so {e[a']\e[a]) = 1 . As a result we get 

P'ic,k) = ^,\ E exp[^]r 

a:x'^=x^mod.{n) 



= P{c, k) 



(2.11) 



For the second case, suppose that the environment 
causes a complete decoherence. If we indexed the ele- 
ments of the reduced density matrix by a and a', then 
this means that its off-diagonal elements vanish. Such 
a case has been formulated in the dynamic theory of 



of the decohering factors F^{a,a') = {elUl} {t)U^it)\e} 
with norms less than unity. In the macroscopic limit 
TV oo, it is possible that F{a, a') 0, for a' ^ a, 
namely, (e [a'] I e [a]) = 6aa'- Then we have 



p'{c,k) 



k)/r] 



<4^ = - 
q'' r qr 



(2.14) 



and 



p'{s) = r(j){r)p' {c, k) 



(j>{r)/q < (i>{r)/n^ <l/n 



(2.15) 



In the remaining part of this section let us proceed on 
to discuss the possible influence of environment on the ef- 
ficiency of Shor's prime factorization algorithm. Generally 
speaking, a deterministic algorithm is said to be efficient 
if the number of the computation steps taken to cxcute 
it increases no faster than a polynomial function of In N 
where N is the input. For a randomized algorithm this 
definition should be modified to fit in the probability 
character. Suppose the one-try-success probability of an 
randomized algorithm ^ is s, then A is said to be efficient 
if Ve > 0, 3 p{x) such that VA^ • (1 - .s)p(1"^) < e, where 
p{x) is a polynomial. Obviously, the polynomial p{x) 
here should have real coefficients and satisfy p(ln A'') > 0. 
All the polynomials appearing in the following are tac- 
itly assumed to have this property. It is also clear that 
in quantum computations all algorithms should be ran- 
domized ones. 

Let A be a quantum algorithm. Suppose for an input 
N the one-try-success probability of A is f{N) where / is 
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a real continuous function defined on the real line. Then 
we have the following lemma. 

Lemma. If there exists a polynomial p{x) such that 
lim (1 - /(iV))Pa"JV) < 1 (2.16) 

AT— >C!0 

then A is efficient. Conversely, if for an arbitrary poly- 
nomial p{x) we have 

lim (1 - /(iV))Pa"JV) > 1 (2.17) 

Af— ►oo 

A is not efficient. 

Proof. Let p{x) be a polynomial such that 
limAr^oo(l - /(iV))P(i'^^) < 1. Then Ve > 0,3a(e) 
such that ( limjv^oo(l - /(iV))p(i"^))«(^) < s. Namely, 
limjv^cx>(l - /(Ar))p(i"^)«(^) < s. Defining 

p'{x) = a{e)p{x) 

we come to the conclusion that there exists some A^o such 

that ViV > No, (1 - /(Ar))p'(i"^) < e. It is now evident 
that one can choose a suitable polynomial ^'(a;) such that 
V7V, (1 - /(A/'))?''"^) < s. This proves the first part of 
the lemma. 

For the second part of the lemma, if the conclusion 
were not true, Ve > 0, there would exist a polynomial 
p^{x) such that ViV, (l-/(7V))P=(i"^) < e.Thus we would 
have limAr^oo(l - /(7V))p=(1"^) < e,leading to the con- 
tradiction 1 < £. The lemma is consequently proved. 

Before concluding this section let us take A to be Shor's 
prime factorization algorithm and return to the above 
mentioned two extreme cases. For the first case, we have 
f{N) > 1/3 In AT. As a result. 



lim (1 - f{N)) 



SlnJV 



< 



lim (1-1/(3 In AT)) 

iV— >cx) 



1/e < 1 



SlnJV 



(2.17) 



So according to the lemma A is efficient. For the 
second case,we have f{N) < 1/N. It is easy to 
prove lim]v^oo(l — l/N)^'^ ^ = 1 for all ntegers m 
so for all polynomials p{x) limAr^oo(l — 1/N)p^^'^^^ = 
1. Consequently, for all polynomials p{x) 



lim (1 

JV->oo 



p{lnN) 



> lim (1- l/iV)P(''^^) = 1 

JV— >oo 



(2.17) 



This means that the algorithm is no longer efficient in 
this case. 



3. Dynamic Model for Decoherence in Quantum 
Computation: 
Generalized Hepp-Coleman Approach 



In this section we begin to study decoherence problems 
caused dynamically by environment. Since the influence 
of environment on quantum computations is reflected in 
the decohering factor (e[a']|e[a]) caused by coupling, as is 
shown in the last section,our discussion will concern the 
microscopic dynamics of the interaction between a qubit 
system and the environment. It is modeled in terms of 
a generalization of the HC model of the wave function 
collapse in quantum measurement. 

In the generalized HC model [30-32 ] , the environment 
E or the detector is made up of N particles and has the 
free Hamiltonian in a general form 



N 



D 



(3.1) 



where the single particle Hamiltonian only depends 
on dynamical variables Xk (such as the canonical coor- 
dinate , momentum and spin ). Let \n){n = 1,2,...L) 
be states of the system (e.g., a quantum register) corre- 
sponding to energy levels £J„ (n = 1, 2, ...M). The system 
S with Hamiltonian 



H, = Y,En\n){n\ 



(3.2) 



n=l 



interacts with the environment through the quantum 
non-demolition (QND)[43] interaction 



Hi = ^^gn,i{xj)\n){n\ 



(3.3) 



n J 



In a special realization for quantum computation, |n) 
may denote an array of qubits, each of which has 
two states |0) and |1), through the definition \n) = 
\non\n2,..., ul-i) where the labels satisfy the unique bi- 
nary representation 



L-l 

n=^n^2\n^=0,l). 

i=0 



(3.4) 



It should be emphasized that, for quantum measurement, 
the interaction between S and E must be chosen to have 
the different strengths for the different states of S, i.e., 
it is required that Qnj ^ Qruj for m ^ n. In fact, the 
so-called measurement is a scheme to read out the states 
of S from the counting number of detector such that dif- 
ferent counting numbers should correspond to different 
states of S. However, this requirement of non-degeneracy 
is not necessary when we extend this generalized HC 
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model such that it is applicable in quantum computa- 
tion. 

If the coupling of the system to the environment is 

degenerate, namely, gn,j{xj) = grn.j{xj) for certain n ^ 
m, we can group the coefficients of the interaction as 
follows 

9i,j = •••• = 9dij = 
Qdi+ij = •••• = 9di+d2j = K2,j, 



5<ii + ..+dg-i + lj — •••• — 9di + ...+dgj = Kqj 

(3.5) 

Correspondingly, the Hilbert space V : {|?^)|'^ = 
1, 2, ...L}of S is decomposed into a direct sum 

q 

of the subspaces 

: {\n = m) = \\,m)\m = 1, ...di} 

: {|n = m + di) = |2,to)|to = 1, ...^2} 



V : {\n = di + .. + dq-i +m) = \q,m)\m = 1, ...,dq} 

(3.6) 

Notice that Kgj 7^ Kq'j for q q'. With the above decom- 
position of the Hilbert space, the interaction Hamiltonian 
can be re- written as 

Kqj {x j)\q,m){q,m\ (3.7) 

g,m j 

We observe that the coupling has the same strength for 
the states belonging to the same subspacc V^. 

In some cases the above classification of state vectors 
is a reflection of the structure of some irreducible rep- 
resentation of certain group chain G Z) K where G, K 
are such chosen that Hi is G— invariant and is at 
most K— invariant. For instance, consider the group 
chain SO{i) D SO (2) . It defines the standard angu- 
lar basis I J, M) through the Casimir operators and 
J3 of 50(3) and 50(2). In this case, we can take 



Hg = H{,P , Js) to be the Zecman Hamiltonian in a cen- 
tral force field (not Coulomb field) if the interaction Hj 

A special cases of the above general discussion has al- 
ready been given in ref. [21]. They introduce a totally 
factorized interaction as of form Hj = Q ^ X^i^i fji^j)- 
Here Q is a system variable commuting with the free 
Hamiltonian Hs of the qubit system and Xj (j=l,2,....N) 
are the variables of the environment with the free Hamil- 
tonian Hjj = X^fe' ^k- The Hilbert space Vg for the sys- 
tem can be spanned by \q, m){m — 1, 2, rf^ for a given 
q ) , the common eigenstates of Q and Hs labeled by q 
and m. 

Q\q.,m) = eq\q,m),Hs\q,m) = Eqm\q,m) (3.8) 

In this case we have the direct sum decomposition Vs = 
J2q®^'^ ^ith the eigen-spaces = Span{\q,m)\m = 
1, 2, dq}.TMs special interaction can be extended to a 
most general form with many system variables Qj{j = 
1,2..K) that cancel certain subspaces of the qubit sys- 
tem simultaneously [21 ]. In the view point of group rep- 
resentation theory, this generalization enjoys an elegant 
mathematical structure [22, 23]. 

4. State Reduction in Time Evolution 

In this section it is shown that the above general struc- 
ture of subspace decomposition (3.5) indeed dynamically 
leads to a scheme grouping the states of qubit system to 
avoid decoherence. Let = Vi ® V2 ® .... ® V/v-i ® Vn 
denote the direct product Hilbert space for the environ- 
ment. Here Vk(k =- 1, 2, denotes the Hilbert space 
of the k'th particle comprising the environment. We will 
prove that, with the QN interaction (3.6) any coherent 
superposition Cm j?, rri) of the states belonging to the 
same subspace is decoherence free while the coher- 
ent superposition ^gk' ™<?)of states belonging to 
difl[erent subspaces may experience a WFC or decoher- 
ence. For a given initial state |cr(0)) = |crj(0)) S Vd, 
of the environment, and a given initial state |/(0)) = 
XmqC^lg, m), of the system , the general initial state 
of the total system |i>(0)) = |/(0)) ®\o{^)) wiU evolve 
into an entangling state 

q m 

^UeM-itlHk + i^gj{xj)]}\aj{0)} (4.1) 
j 

where Eqm = Edi+..+d,_i+ni,rn = l,2,..dq. Notice that 
the states belonging to a subspace entangle with the 
environment through the same factorized components 
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l[Ujit)\ajm ^ 

3 

nexp{-*[i7fe + K,,j(x,)]t}|a,(0)) (4.2) 



This kind of factorization structure in evolution of wave 
function is crucial to decoherence or WFC [30-36]. In 
fact, the reduced density matrix of the system at time t 
is 

p{t)=Tr,mt))m)\) = 

Y,{Y.\CU''\q,m){qM + 

q m 

^ exp[iEgm,t - iEg^tjC^C!^, \q, m){q, m'\} 

q^q' m.m' 

N 

{q',m'\ X{{aAm]'\t)U]{t)\aM) (4-3) 

where Tr^ means taking partial trace over the variables 
of the environment. Prom this expression we see that 
each off-diagonal element of , labeled by q and q' , is 
accompanied by a decohering factor 

JV 

in the form of factorized function .Obviously, if the 
initial state |/(0)) belongs to a single subspace V^, 
then the terms accompanied by F^^^i {N, t) do not ap- 
pear. Thus the system will remain in the pure state 
eyip[-iHst]\f{0)){f{0)\exp[iHst] throughout the evolu- 
tion process. This fact is significant for developing 
schemes to carry out error free quantum computations. 
The expression also manifests the happening of deco- 
herence for those states belonging to different subspaces 

Next we consider the dynamic process of decoherence 
when a superposition of states mixes the vectors belong- 
ing to different subspaces. Naively, as F^j^^>{N,t) is a 
multiplication of N factors ^ it) with norms not larger 
than the unity, it may approach zero in the macroscopic 



limit with very large A^. To deal with this problem pre- 
cisely, we define an real number not less than zero 

Af'(i)=-ln|F;;^,(t)| (4.5) 

Then the norm of the accompanying factor F^^^ {N, t) is 
expressed as 

N 

\F^^g>{N,t)\=eM-Y.^T'\ (4-6) 

Obviously, the series Yl^=i (*) — since each term 
is not less than zero. There are two cases in which the 
accompanying factor F^^^' {N, t) approaches zero in the 
macroscopic limit with very large N. The first case is 

that the series J^JLi^f^ (*) diverges on (0,oo]. The 

second case is that the series converges to a monotonic 
function which approaches zero as iV ^ oo. Therefore, it 
is possible that 

q m 

exp[iEgm't-iEgmt]ClCZ\q,m){q,m'\} (4.7) 

as N oo. In the next section, some examples will 
be presented to illustrate the above mentioned circum- 
stances explicitly. 

5. Dynamic Decoherence in Environment 
Consisting of N Two- Level Subsystems 

To make a deeper elucidation of the above mentioned 
quantum dynamic mechanism of decoherence in quan- 
tum computation and the relevant scheme of grouping 
the states of a qubit system to be prevented from de- 
coherence, in this section, we model the environment as 
consisting of N two level subsystems. We recall that 
Caldeira and Leggett [37] have pointed out that any en- 
vironment weakly coupling to system may be approxi- 
mated as a bath of oscillators . On the condition that 
"any one environmental degree of freedom is only weakly 
perturbed by its interaction with the system" , they have 
also justified describing the influence of environment by a 
coupling linear in the bath variables up to the first order 
perturbation. With this justification, we observed that 
any linear coupling only involves the transitions between 
the lowest two levels (ground state and the first exci- 
tation state) of each harmonic oscillator in environment 
though it has many energy levels. Therefore in this case 
we can also describe the environment as a combination 
of many two level subsystems without loosing general- 
ity. In fact, for quantum computation, Unruh[ll] and 
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Palma et al [12] have considered the harmonic osciUa- 
tor environment. Their model is equivalent to one in- 
troduced for the WFC in quantum measurement by Sun 
et.al [32, 31,3 4] . A similar model has also been touched 
by Leggett et al [37,38] and Gardiner [39] in studying the 
tunneling effect in a quantum dissipative process. Here 
we choose cquivalcntly the two level subsystem model to 
manifest some characters independent of environment in 
the weakly coupling limit and to demonstrate explicitly 
the qualitative calculation of decoherence time through 
a sample example without quantum dissipation. 

Let \gj) and \ej) be the ground and excited states of j 
'th subsystem .We define the quasi-spin operators 

= |ej)tel + \9i){ej\ 

CT2(j) = -i[\ej){gj\ - \gi){ej\] 

(^3{j) = \ej){ej\ - \gj){9j\ (5.1) 
Then we introduce the Hamiltonian of the environment 

N 

He = J2ncOja3{j) (5.2) 
and the interaction coupling to a qubit system 

JV 

Hj = f{S)J2^gja2{j) (5.3) 

where f{S) is function of the variable S of the qubit sys- 
tem. 

In this section let us mainly focus on the simplest case 
where the system consists of two qubits with the Hamil- 
tonian 

Hs = hviSsil) + n7j2S3{2) (5.4) 

where S{1) = (8) 1, Ssi^) = 1 (g) cr^, (s = 1, 2, 3) denote 
spin operators acting on the first and the second qubits 
respectively; (Js(s = 1, 2, 3) denoting the usual Pauli ma- 
trix. We consider the special interaction given by 

/(5) = 53(1) + 53(2) (5.5) 

It means that in our model the interaction has the same 
strength for different states. This model is very simple, 
or even too simple in some sense. But we would like to 
point out that the so called Free Hamiltonian Elimination 
model in ref.[21] is substantially only a plain generaliza- 
tion of the present example to the multi-pair case if one 
takes into account the SU (2) rotation transformation. 

Let |1) and |0) be the qubit states that satisfy Sslk) = 
(— )'^+^|fc), (fc — 1,0). With the chosen interaction form, 
the Hilbert space, spanned by 



{|1,1) = |1)®|1),|1,0) = |1)®|0), 
|0,1) = |0)®|1),|0,0) = |0)®|0)} 

contains a null subspacc V'^ of Hj spanned by 1 1 , 0) and 
|0, 1) . Any superposition |0(O)) = A\1,0) + B\0,1) in 
this subspace will preserve its purity in evolution pro- 
cess though the system has interaction with the envi- 
ronment. Precisely, the pure state 10(0)) ((/)(0)| will 
evolve into the pure state Uo{t)\(j){0)) {4>{0)\UQ{t) where 
Uo(t) = cxp[— 277ii53(l) — 2772^53(2)] is the free evolution 
operator of the qubit system. Physically, this fact implies 
that no useful information leaks out of the system in the 
process and the coherence is preserved. This analysis can 
be easily generalized to the many bit case where the free 
qubit Hamiltonian takes the form = J2k=o ^VkSsik) 
and its interaction with the environment is determined 

by 

L-l 

/(5) = ^A,53(fc)] (5.6) 

fe=0 

where L is the number of qubits 

k—ltimes 

used and Ss{k) =1 (g) (g)l (g)(7s (g) 1 (g) • • • » l.The dif- 
ferent A^s indicate that each single qubit has a different 
coupling to the same environment. In the Hilbert space 
of this L-qubit system with the basis 

|a> = |ao> (g ||gi) (g jga) (g) ... (g) jgi-i), 

gfc = 0,1; fc = 0,1,2,. ..,L (5.7) 

the subspace preserving coherence can be spanned by 
those basis vectors satisfying 

L-l 

J2 = const.e (5.8) 

fe=0 

Let us return to the two qubit example. If a super- 
position contains a vector outside the decoherence free 
subspacc, decoherence will happen in an entanglement of 
system state with environment state. For example, if the 
initial state \ip{0)) = C|0,0) + L»|l, 1) of the system in- 
volves states not belonging to while the environment 
is initially in the vacuum state |0)e = |fli)(g'|52)®""® j^jv) 
where \gj)is the ground state of j'th two level subsys- 
tem , the corresponding pure state density matrix |</3(0)) 
((p(0)|(g|0)e(g)e(0| of the total system formed by the qubits 
plus the environment will experience a unitary evolution 
to reach a pure state prit)- Its reduced density matrix 

p{t) = TrePTit) = \C\^\0, 0)(0, 0| + \D\^\1, 1)(1, 1| 
+[CD* exp[2i(r?i + r]2)t]F{N,t)\0, 0)(1, 1| + H.c] (5.9) 
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is no longer pure because the environment state becomes 
correlated with the system state. Here the decohering 
factor 

N N 

F{N,t) EE n F,{t) = l[{9j\Uj,mjo{t)\gj) (5.10) 
is determined by the effective evolution operators 

Uja{t) = exp[-iWj(T3(j)t - i^agj(r2{j)t] 



a =2,^0 = -2, (a = 0,1) 



(5.11) 



corresponding to the qubit states |0, 0) and |1, 1) respec- 
tively. Using the formula exp[ia^ ■ A] = cos A + ilf ■ 
nX sin A for a given vector A of norm A along the direc- 
tion nX, we get the explicit form of Uja{t) 

Uja = COs{fljat) — i [a2{j)sm9ja + a3{j) cos 9 ja]sm{Q,jat) 

(5.12) 

where tan 9ja = %^ , ^ja = -^{gjCay + ^j- Then, we 

JV 

get the decohering factor F{N,t) = fl Pih't), which 

i=i 

is an A^— multiple product of the factors F{i. f) = 1 — 
2 sin^ 9j s\v? Qjt of norm less than 1. Here, we have used 



the definitions tan 9j 



2s, 



Agj"^ + LuHoT the 



special labels ^ = 2,^o = —2. Therefore, the temporal 
behavior of decoherence is described by 



TV 



\F{N,t)\ =exp^ln|l-8gsin2(17jt)| (5.13) 



i=i 



In the weakly coupling limit that gj <C u>j, we get 



N 



\F{N,t)\ = e-^W = exp - ^ ^ sin^ {u;,t) (5.14) 



A special case is that the subsystems constituting 
the environment are identical and the environment has 
a constant discrete spectrum, i.e., ojk =constant oj, 
gk =constant g. In this case, the off-diagonal elements 
with the factor exp{— ^^|- sin^ wtjapproach zero as ^ 
oo for all t except those satisfying wt = 2k'K{k = 0, 1, 2...). 
For general information, one needs a detailed analysis 
about the behavior of the series S{t) for various spec- 
trum distributions of the environment. Of special interest 
is the case with continuous spectrum. In such case S{t) 
can be re- expressed in terms of a spectrum distribution 
p{ijJk) as 



S{t) 



~ Jo <4 



P{'^k)gk sin^ uJkdoJk 



(5.15) 



Notice that, in the case of discrete spectrum, the dis- 
tribution means a degeneracy : there are p(w/c) subsys- 
tems possessing the same frequency u)k- Prom some con- 
crete spectrum distributions, interesting circumstances 
may arise. For instance, when p{ujk) = -^l/gk the in- 
tegral converges to a negative number proportional to 
time t , namely, S{t) = —"/t. This shows that the norm 
of the decoherence factor is exponentially decaying and 
as f — > oo, the off-diagonal elements of the density ma- 
trix vanish simultaneously! Another example of contin- 
uous spectrum is the Ohmic type [37,38] p{u)k) ~ — 
which leads to a diverging integral S{t) oo for t 0. 
In conclusion, in the present example, we can choose a 
suitable spectrum distribution of the oscillators in the 
detector, such that the series S{t) diverges to infinity, or 
in other words, the dynamical evolution of the system 
plus environment results in the complete decoherence in 
the reduced density matrix of S independent of the tem- 
perature. However, it is only an accident situation owing 
to the special choice of the initial state. For a general 
initial state, we will see, the decoherence process indeed 
shows a temperature independence. 



6. Decoherence Time for L-Qubit System 



Usually, when the coherence of a quantum system de- 
velops a characteristic decay proportional to a factor of 
the form exp(— t/t^), ta, which characterizes the speed 
of the decoherence or the transition of the system from 
the quantum regime to the classical one, is called the 
dec;oher(;nc(; time. Its value depends on the physical fea- 
ture of the quantum system and their interaction with 
the environment. For a single qubit system some numer- 
ical estimates of td have been made by DiVinccnzo[44] 
for several physical realizations. It ranges from lO** s (for 
nuclear spins )to 10~^^ s (for the electron-hole excitation 
in bulk of a semiconductor). In practice, to carry out 
a quantum computation, one needs a large number of 
qubits, e.g., in Shor's algorithm factoring large number 
n, i oc Inn. Accordingly, in the following we extend the 
dynamic analysis to show how the speed of decoherence 
becomes larger as the number of qubits increases. 

Let us consider the L-bit system coupling to the envi- 
ronment mentioned in the last section. The interaction 
constants A/j's are chosen so that eigen- values 



L-l 



(6.1) 



fc=0 



are not degenerate for {q^ = 0, 1}. Starting from an initial 
state \ip{0)) = C\q) + D\q'), where 



L-l 



\p) = Jl <»\Pk),p = q,q' 



(6.2), 



k=l 
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the initial pure state density matrix of the total system 
formed by the qubits plus the environment will experi- 
ence a unitary evolution to reach a pure state pL{t)- 

Imitating the calculation process in the last section, wc 
can obtain the reduced density matrix p{t) = TrePL{t)- 
Its off-diagonal elements are proportional to the decoher- 
ing factor 

N 



Consequently 



FL{N,t) = l[FL{j,t) 



N 



]l{gj\UliL,t)U,,,{L,t)\gj) (6.3) 



where 



N 



Ujq{L, i) = n exp[-ia;jO-3(j> - ii{q)gja2{j)t] (6.4) 

Using the notions 

tane,(g) = ^, %(5) = Jfe^(g)P+a;| (6.5) 

and the matrix representation of Ujq{L, t), after straight 
calculation we get 

Fl {j, t) = sin {q) sm[Q.j {q)t] sin {q') sm[Q.j {q')t] 



+{cos\p.j{q)t] — icos9j{q) sm[Qj{q)t]} 



X {cos[n j{q')t] +icosej{q') sin[%(g')i]} (6.6) 

Trivially, F]^{N,t) becomes unity when q = q' ■ How- 
ever, when q ^ q', in the weakly- coupling limit gj <C u)j, 
we have 

sm0j{q) ~ ej{q) 



cosej{q):^l--e]{q) 



flj (q) ~ LOj 



(6.7) 



Thus FL{j,t) ~ 1 - ^{9,{q) - ej{q')}^sm\ojjt) + 
" dj{q')} sm{2ijjjt). Since in such weak-coupling 
limit Oj{q) ~ sin 6 j{q) ~ ^{q)^, we obtain the decoher- 
ing factors 



FL{j,t) 1 - |^U(«) - UllV^'^^'i^jt) 



• 2 

+S{^'(9)-^|(5')}sin(2a;,0 



(6.8) 



|Fz,(iV,i)|=exp{-[e(g)-e(g')] 



N 2 

2 \ " yj 



sin^(ci;ji)} 
(6.9) 

In summary, the temporal behavior of the dccoher- 
encc is described by F{N,t), and actually determined 
by \F{N,t)\ , which is of the form exp[— S'L(i)]- Here 

SL{t) = m) - ai')? Ef=i 4^ sm^u;,t). For identical 
qubits Afe =l,the fastest decoherence happens between 
the two initial states \q) = |go = 1) 'X' |gi = 1) (E" • • • ® 
\qL-i = 1) and \q') = \qa = 0)®|gi = 0)(^---® = 0) 
. In this case \FLiN,t)\ = exp[-L^S{t)]. Thus for the 
instance with S{t) = jt, which is discussed in the last 
section, we have \F{N,t)\ = cxp[— L^7i] where 7~-^is the 
decoherence time for a single qubit. This shows that the 
cliarac;tc!rized time of the fastest decoherence happening 
in the L-qubit system is times of that of a single qubit. 
This conclusion first obtained by Palma et.al.[12] is given 
here in the framework of quantum dynamic model of de- 
coherence. 



7. Temperature Dependence of Decoherence 



The above discussion about decoherence in quantum 

computation only concerns the situation of zero temper- 
ature. In this section we consider the influence of en- 
vironment at a finite temperature. Suppose the initial 
state of the total system is described by a density ma- 
trixp(O) = Ps(0)®pb(0) where p,{0) =|</)(0))((/.(0)| is the 
density matrix of the system, while Pb{0) is that of the 
bath 



Pb{0) 



Trbexp{-(3Hb) 



N 



N 

n 2 cosh (Pcoj) 



(7.1) 



with /3 = 1/{KbT). For the initial state |0(O)) = 
A|0,0) + we obtain the same decohering fac- 

tor F2{N,t) = nf [1 - 2sm^ ejsin'^injt)] which is in- 
dependent of temperature. This is due to the special 
choice of the initial state with a certain permutation 
symmetry between |1,0) and |0, 1). For a general initial 
state |i^(0)) = C\q) + D\q'), we can calculate the factor 

FLiN,t) = n FL{j,t) as follows 



FLij,t) 



Trt,[Uj,,{L,t)pjbmUL,t)] 



' sin 9j (q) sin[Oj {q)t] sin 9j {q') sm[Qj {q')t]- 
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cos[Slj{q)t\ cos[Clj{q')t\ + 
cos 9j (q) sin[f2j {q)t] cos 9j {q') sm{Qj {q')t) 

-|tanh(/3wj)sin(2f7j((7)t){cos0j((7) - cos0j{q')} (7.2) 

Notice that the effect of finite temperature only ap- 
pears in the imaginary part of the decohering factor. 
In the weakly coupling limit, it is not difficult to ob- 
serve that FlU, t) ~ 1 - ^{Oj{q) - Oj{q')y sm^[u)jt] - 
lt£inh{Pujj)sm[2ujjt][9]iq) - 6'|(g')] or 

FlU, i) ^ 1 - ^{^(9) - ^j{q')r sin2(a;,-i) 

■ 2 

- ^1(9')} tanh(/3a;,) sin(2a;,t) (7.3) 

It reflects the novel fact that, in an environment weakly 

interacting with the qubit system, the decoherence time 
do not depends on temperature as a result of the 
temperature-independent norm of FL{j,t). Thus in this 
case thermal fluctuation plays a role in quantum com- 
putation only through affecting the phases of the off- 
diagonal elements of the reduced density matrix. 

8. Universality of Environment in Weak Cou- 
pling Limit 

An environment surrounding a qubit system for quan- 
tum computation maybe very complicated. Intuitively, 
the dynamic process of decoherence in quantum compu- 
tation should depend on the details of interaction be- 
tween the qubit system and the environment. So gener- 
ally it seems impossible to control decoherence in a qubit 
system. Nevertheless, one may well expect that in some 
limit situations there exists certain universality in the 
dynamics of interaction so that the physical parameters 
(such as the decoherence time and decoherence factor) 
dominating a quantum computation process would not 
depend on the detail of environment . For the tunnel- 
ing problem in quantum dissipation process, this kind of 
universality has been considered by Caldeira and Leggett 
[37, 38] by modeling the environment as a bath of har- 
monic oscillators with a linear coupling to the system. 
In this section, we illustrate that, in the weakly coupling 
limit, the above results obtained from the two- level sub- 
system model of environment coincide with those from 
the harmonic oscillator model concerned in various quan- 
tum irreversible processes, such as wave function collapse 
[32,34,31] and quantum dissipation [37-42]. 

Let a^and Oj be the creation and annihilation op- 
erators for the i'th harmonic oscillator in the environ- 
ment. The Hamiltonian of the environment takes the 



form H = X]j=i iT-'^jCij'ci and its interaction with the 
qubit system can be modeled as a linear coupling: 

JV 

Hj = f{s)J2^gj{af+aj) (8.1) 

where f{s) is a linear or non- linear function of the qubit 
system variable s. Let the initial state of the qubit sys- 
tem |</j(0)) = A\a) + B\P) be a coherent superposition 
of two eigenstates of s,s\a.) = a\a),s\i3) = f3\f3) and let 
the environment be initially in the vacuum state |0)e = 
|0i) IO2) (8) • • • \0n) where \Oj)is the ground state of 
the i'th single harmonic oscillator. The corresponding de- 

cohering factor F(iV,t) = nf=i h{0\U-\t)U^{t)\0)h = 
Ylj'=iFj{t)c&u be obtained by solving the Schrodinger 
equations of Uj{t) (7 = a. (3) governed by the Hamilto- 
nian of forced harmonic oscillator 

Hjt = nujjofaj + f{j)gj (a+ + aj) (8.2) 

In fact, by the so c;alled Wei- Norman algebraic expansion 
technique one has the following explicit result [12,32,35]. 

F{N, t) = exp{-[/(a) - f{(3)f M siu^"^)} 
X eM-maf - m'] E ^[t + ^^]} (8.3) 

— LVj ivj 

J = l ■' ■' 

The decoherence time is decided by the real part of 
F(N,t), which is the same as in eq.(6.8) from the two 
level model of environment in the weakly coupling limit. 
This can easily be seen if only one replaces ^ in the 
above equation by LOj . Thus in the weakly coupling limit 
the differences among different models of environment 
arc only reflected in the imaginary parts of the decoher- 
ing factors This simply implies that in this case the de- 
tails of environment does not affect the speed at which a 
quantum system approaches the classical kingd(mi. But 
they do affect the success probability of a quantum com- 
putation. 

9. A Discussion 

We have seen from our model that for a quantum regis- 
ter with L qubits the relevant coherence develops a char- 
acteristic decay proportional to cxp[— L^/f^f] where t,i is 
the typical decoherence time for a single bit. Thus if a 
quantum algorithm calls for K elementary computation 
steps and each step takes time r on the average in order 
that the algorithm could be feasible we should have the 
condition 

L\K < td (9.1) 
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Generally speaking, this would pose a strong restriction 
on L and K. We need to develop proper quantum er- 
ror correction schemes to cope with this difficulty caused 
by dccohcrcncc , which is unavoidable in the quantum 
kingdom. Along this line there has been some progress. 
Nevertheless there is another severe problem which may 
endanger the assumed great utility of quantum comput- 
ers. In section 2 we have shown that environment may 
affect the efficiency of a quantum algorithm. Although 
our discussion is not sophisticated enough it indeed gives 
us a frustrating information. This problem deeply rooted 
in the quantum kingdom seems to have been ignored. We 
think it is now time to face it seriously. 

This work is supported in part by special project of the 
NSF of China and the National Excellence Youth Foun- 
dation of China 
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